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ABSTRACT 

We present 3" resolution imaging of the z = 4.7 QSO BR1202-0725 at 900 /im from the Submil- 
limeter Array. The two submilhmeter continuum components are clearly resolved from each other, and 
the positions are consistent with previous lower frequency images. In addition, we detect [C II] line 
emission from the northern component. The ratio of [C II] to far-infrared luminosity is 0.04% for the 
northern component, and an upper limit of < 0.03% is obtained for the southern component. These 
ratios are similar to the low values found in local ultraluminous galaxies, indicating that the excitation 
conditions are different from those found in local field galaxies. X-ray emission is detected by Chandra 
from the southern component at Lo.5_2feey = 3 x lO"*^ erg s~*, and detected at 99.6% confidence from 
the northern component at Lo.5_2fcey ~3xl0'''*erg s~*, supporting the idea that BR1202-0725 is a 
pair of interacting galaxies at z = 4.7 that each harbor an active nucleus. 

Subject headings: galaxies: formation, galaxies: starburst, cosmology: observations, galaxies: high 
redshift, submilhmeter, galaxies: individual (BR1202-0725) 
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1. INTRODUCTION 

The far infrared (FIR) fine structure lines of abun- 
dant species (e.g. [C I], [C II], [O I]) are recognized as 
one of the most important coolants in the warm ISM. In 
particular, the emission from singly ionized carbon (i.e. 
[C II], Xrest = 158/im) is known to trace warm (> 200 



K) and dense (rir 



jhoto-dissociation 



regions (PDRs; iKaufman eranil999|) . directly identi- 
fying the gas surrounding active star forming regions. 
The [C II] line is the dominant coolant in the Galac- 
tic ISM, accounting for 0.1 - 1 % of the FIR luminosity 
l|Hollenbach fc Tielen3ll997f) . The brightness of the line 
makes it an attractive tracer of star formation in galax- 
ies, and a relation holds between [C II] luminosity and 
the massive star formation rate derived from Ha e mission 
for ga laxies that have Lfir < 3 x 10*° (Boscli r"et al.l 
I200a) . In the FIR luminous galaxies (ULIRGs), the 
[C Il]-to-FIR flux ratio is an order of ma gnitude smaller 
than in normal star forming galaxies (Luhma n et al.l 
[1^98, 2003). The reason for the deficiency in the 
[C II] line luminosity is not well unders tood, but sev- 
eral possible scenarios have been offered (iLuhman et al.l 



erai p ossible scenarios nave been onerea (iMmiliaJi et_au 
19981 IMalhotra et alJl200lt IPierini. Leech fc Volkll200a 



Luhman et al]l2003|) . 



Despite the large line luminosity expected, detect- 
ing [C II] line emission from luminous high redshift 
sources using a single- dish submil l imeter telescope has 
proved difficult (e.g. Ilsaak et alJ Il994t Ivan der WerU 

m 



lltMlBolatto et al.ll2004HMarsden et al jl2OO50. The 1^ 

[C II]-to-FIR ratio suggests that the physical properties 
in high-z sou rces resemble thos e of the local ULIRG pop- 
ulation fe.g. iBlain et alJl2002() . Detection of the [C II] 
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fine in the z = 6.4 QSO J1148+5251 has been reported 
recently bv lMaiolino et al . ( 2005). Here we report a first 
detection of the [C III line emission from the z — 4.7 
hyper-luminous (- 10" L©) QSO BR1 202-0725 observed 
the Submillimeter Array (SMA; iHo. Moran fc Lol 

5 



BR1202-0725 is an optically bright QSO {Mb = 
28.5) with a 2" 3 Lyman a extension t o ward the north- 
west llHu. McMahon fc Egamil 119961 iPetitiean et al 



19961 iFontana et al.lll998t lOhvama. Taniguchi fc Shiova 



2004[) . Radio continuum emission was detected at 1.4, 
4.9 and 43 GHz (Yun ct al. 2000; Carilh ct al. 2002). 
The emission at 1.4 GHz was resolved into two com- 
ponents; one at the location of the QSO and the other 
^ 4" northwest. This apparent co mpanion is also seen 
in millimeter continuum emission fOmo nt et all Il996t 
iGuillotcau ct al. 1999). Abundant CO emission (Mjy^ > 
10*° Mq) was detected toward bot h the opt ical QSO 
and the northw estern companion fOmont et al.. ■199§ 



Carilli et al 



.1 120021) 



lOhta et alJ IT99'6 : GuiU oteaul l200lt 
The CO-derived rcdshifts are z = 4.6947 and 4.6916 for 
the QSO and the northwestern companion, respectively 
ijOmont et al.llT99Q) . 

We adopt = 70 km s"* Mpc"*, ri,„ = 0.3, Qa = 
0.7, which gives Dl = 43.27 Gpc (1" = 6.5 kpc) for the 
luminosity distance to BR1202-0725. 

2. OBSERVATION AND DATA REDUCTION 

A total of 5 tracks (Table ^ were obtained toward 
BR1202-0725, where the total on-source integration time 
was 19.6 hours. The SIS receivers were tuned to a red- 
shifted [C II] frequency of 333.969 GHz (900 /im) (LSB). 
The SMA correlator had 2 GHz total bandwidth with 
0.8 MHz (0.7 km s~*) spectral resolution. The data 
were calibrated using the IDL based SMA calibration 
tool MIR. The quasars 3C279 and 3C273 were observed 

^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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to calibrate the antenna based time dependent gain, and 
Callisto, 3C279, 3C273, and 1924-292 were used to cali- 
brate the bandpass. The flux scaling is derived assuming 
Callisto has a total flux density of 15 Jy at 334 GHz, and 
we estimate a 20% uncertainty in the derived fluxes. 

Imaging and analysis of the visibility data were carried 
out using MIRIAD. The angular resolution achieved by 
adopting ROBUST = 2 was 3'.'4 x 2'.'7 (P. A. = -8°). The 
continuum image constructed using the line-free data has 
an RMS noise of 3.0 mJy. Although the varying S/N 
across the the SMA bandpass makes it difflcult to char- 
acterize the overall noise properties, the RMS noise in the 
[C II] spectrum calculated using several line-free regions 
of the data cube is estimated to be 7.5 mJy. 

3. RESULTS 

Figure ^ (left) presents the 900 /im continuum im- 
age of BR1202-0725. The coordinates of the two com- 
ponents are a (J2000) = 12''05"22.98^ and S (J2000) = 
-7°42'30'.'0 for BR1202N, and a (J2000) = 12''05'"23.12'* 
and S (J2000) = -7°42'32'.'5 for BR1202S with uncer- 
tainties of ±0'.'2. The peak of the northern compo- 
nent (B R1202N) is coincident with the CO (2-1) peak 
ijCarilli et al. ,200^1 within the 0'.'2 astrometric accuracy. 
The peak of the southern component (BR1202S) is dis- 
placed to the north from the CO (2-1) and the radio 
continuum peaks (CariUi ct al. 2002) by ~ 0'.'5. The lo- 
cation of BR1202S is consistent with the optical/NIR 
QSO reported by Hu. McMahon fc Egamil l(l9961 The 
derived flux densities from the line-free region of the spec- 
trum are Sgoo = 27 ± 4 mJy (BR1202N) and 32 ± 4 mJy 
(BR1202S). BR1202S emits 54 ± 9% of the total 900^m 
flux, which is slightly lower than the earlier observations 
at 1.3 mm in which 65 ± 10 % of the total flux was found 
in BR1202S ijOmont et alJ [T9961. 

A spectrum at the peak of the 900/im continuum in 
BR1202N is shown in Figure^ (middle). A similar spec- 
trum toward BR1202S is also shown in Figure^ (right). 
Continuum subtraction for both spectra was performed 
by fitting a zeroth order baseline in the line-free chan- 
nels of the visibilities. The BR1202N spectrum shows 
a 22 mJy peak at z ^ 4.691 with a line profile that 
spans the redshift range oi z = 4.687 - 4.696. No 
significant detection is seen in the BR1202S spectrum. 
The FWHM of the [C II] fine is 240 ±50 km s'^, 
which is narrower than the CO ( 5-4) profile ( FWHM 
= 350 ± 60 km s'^) obtained bv lOmont et"all l|l99^ . 
The formal significance of the [C II] line detection over 
the entire line-width shown in Figure ^ is 4(t, assum- 
ing Gaussian statistics. The true significance of this fea- 
ture is likely to be higher because noise is non-uniform 
across the bandpass. The excellent agreement between 
the [C II] line profile and the published CO (5-4) line 
profile further supports the robustness of our [CII] line 
detection. By integrating the emission in the velocity 
range -230 to 130 km s~^, we estimate the [C II] line flux 
of S[cu] — (6.8 ± 1.1) Jy km s^^ which translates to a 
[C II] luminosity of Lpn] = (4.5 ±0.7) x IO^Lq. Using a 
FIR luminosity of Lfir = (1.2 ±0.2) x IO^Lq (see §4.2), 
we obtain L[cii]/Lfir = (3.8 ± 0.9) x lO""* (0.04%) in 
BR1202N. In addition, assuming Av = 500 km s~^ and 
Lfir = 2.6 x lO^'^L© (see §4.2), we obtain a 3a upper 
limit of L[cii]/Lfir < 2.8 x lO^"* (0.03%) in BR1202S. 



4. DISCUSSION 

4.1. Properties of the [C II] Line Emission 

The exceedingly low L[cn]/LFiR ratios for both 
BR1202N and BR1202S are consistent with the simi- 
larly low ratios found in local ULIRGs where the aver- 
age was fou nd to be Ljchj/Lfir = 0.04% (iLuhman et alJ 
ri99Sl I2OO3I) . This is in stark contrast to the values 
found in the local galaxies bv iMalhotra et aTl l)2001t) in 
which ~ 70% of the sources show L[cii]/Lfir > 0.2%, 
with the majority in the range 0.1 — 1%. Past stud- 
ies have reported L[ciii/Lfi r ~ 0.02% for the z — 
6.42 QSO SDSS J1148+5251 lIMaiolino et alJl200l . and 
L[ciii/Lfir < 0-4% i n a z = 4.926 QSO CL 1358+62 
(iMarsden et al.jl20(il . Our new measurement suggests 



^irsQen ' 

that the [C II]-to-FIR ratio (0.04%) in the high-z hyper- 
luminous source BR1202-0725 is comparable to local 
ULIRGs. Table 13 summarizes these results. 

The integrated [C II] line intensity is known 
to be a tracer of star formation activity in nor- 
mal galaxies. We derive the star formation rate 
(SFR) of 2900 Up. yr-i using the calibration of 
iMaio hno etal\ l|2005() . This is 1 - 2 orders of magni- 
tude higher than the SFRs derived from optical recom- 
bination lines observed near BR1202N ( i.e. 10 - 230 
Mp, y r ~^:lHu. McMahon fc Egami 1996: .Petiti ean et all 
19961 lOhta et alJ 120001: lOhvama. Taniguchi fc Shiovi 



2004|) , but consistent with the SFR derived from the SED 



flt to the millimeter/submillimeter continuum measure- 
ments (i.e. SFR ^ 2000 Mq yr'^ in BR1202N; see §4.2). 

Another important diagnostic of the physical proper- 
ties is the [C II]-CO flux ratio. Early studies suggest 
the close spatial and kinematical association between the 
[C II] and CO (1-0) emissions in nearby galaxies, sug- 
gesting that this ratio is a good tracer of star forma- 
tion a£Hvit^;_tlmMsi^ filling fac- 
tor l|Crawford et alJll985l:rStacev et alJll99lD . Assuming 
the luminosity ratio between the CO (5-4) and CO (1- 
0) li nes in BR1202N is similar to that of the Clover- 
leaf (jBarvainis et all IT991 IT991 . the CO (1-0) lumi- 
nosity is estimated to be Lco(i-o) ~ 10^ Lq, yield- 
ing L[cii]/Lco(i-o) ~ 4500. This value is highly un- 
certain since the CO (5-4) - CO (1-0) ratio may be 
different in BR1202-0725. Nevertheless, this is a fac- 
tor of 2.5 larger than the average value s of ULIRGs (i.e. 
L[cii]/Lco(i-o) ~ 1700: Ivan der WerjlT998i). but com- 
parable to the average values of normal galaxies (i.e. 
Lrciii/Lco(i-o) = 4400 - 6300; ICrawford et all IT985t 
iStacev et al.lll99il) . The higher observed ratio may sug- 
gest a factor of few higher Go/n ratio (where Go and n 
are the far UV flux and the cloud density respectively) in 
B R1202N compared to the local ULIRGs (see Figure 9 
of lKaufman et alJlT999(l . 

4.2. Properties of the Continuum Emission 

The spatially resolved submillimeter images allow us to 
investigate the properties of the individual components 
in BR1202-0725. The radio to FIR SED of BR1202-0725 
is shown in Figure |21 where the individual flux densities 
are given for BR1202N and B R1202S for the 1.4 GHz 
fCarim ct al. 2002), 1.3 mm l|Omont et alJ Il996"l and 
900/im observations (sec figure caption for the adopted 
parameters). The SEDs of the two sources are ana- 



lyzed using the model bv lYun fc Carillil l|2002|) . yielding 
Lfir = 1.2 X 1013 and 2.6 x lO^^ ^nd SFR = 2000 
and 4500 Mq yr"! for BR1202N and BR1202S respec- 
tivel y. The discrepa,ncy b etween the SFR derived by us 
and 'Yun fc Carillil (j2Qp)(2300 Mg yr"! for the entire 
BR1202-0725) is due to the difference in the adopted 
cosmology. Note that while the millimeter/ submillimeter 
continuum flux densities of BR1202N and BR1202S are 
comparable (§3), the radio continuum flux density of 
BR1202N is 3 times larger. This suggests a presence 
of a highly obscured AGN in BR1202N (see §4.3), or 
a strong effect of a highly energetic jet from BR1202S 
l|Klamer et al■ll2004^ . 

Assuming a dust a bsorption coefficient of Kd = 
1.5 m ^ kg-i at 160^m fHildebrand' 19831 iDraine k Led 
Il984|) , the dust mass of the individual components is es- 
timated to be Md (BR1202N) = (5.4 ±0.8) x 10^ Mq and 
Md (BR1202S) = (9.2 ± 1.1) x 10** Mq. Using the Hz 
mass derived from the CO (5-4) emission ijOmont et al.l 
Il99(j) . we find MnJMd 50 ± 10 for both compo- 
nents. This ratio is an order of magnitude lower than the 
mean value for LIRGs/ULIRGs (MnjMd = 540 ± 290; 
iSanders. Scoville fc Soifejll99H) . ne arW spiral galaxies^ 
(MHo /Md ~ 500; see Fig. 2b of iDevereux fc Yound 
Il990|) and a sample of Spitze r 160 /xm selected galaxies 
(MnJMd = 330; Yim et aLl l2006D . These gas-to-dust 
ratios are highly uncertain because dust emissivity, tem- 
perature, and CO-to-H2 conversion factor are all poorly 
constrained. If the low ratio for BR1202-0725 is correct, 
however, it may indicate that the H2 mass is significantly 
underestimated in BR1202-0725, and/or the warm com- 
ponent of dust dominates the far infrared emission, sug- 
gesting that the entire ISM is involved in the starburst 
activity. 

4.3. BR1202-0725: A Colliding Galaxy System at 
z = 4.7 

Our analysis of previously unpublished archival Chan- 
dra data (ObsID 3025) clearly shows that BR1202S 
is an X-ray source and gives 26.4 net counts in 
9.64 ksec. The derived flux from BR1202S is /o.5-2fcey = 
6 X 10^1^ erg cm^^ s^^ which yields a luminosity 
of Lo.5_2feey = 3 X 10^^ erg s^^ (assuming F = 
2.0 and Galactic foreground Njj — 3.35 x 10^°cm~^; 
iStark et alJ[T99 ^. X -ray emission is also seen within 
1" of the position of BR1202N. With just 2 counts 

^ The CO to H2 conversion factor used bylOmont et al.' J19961), 
IHandcrs. Scoville & Soifcr 1 1991) and Dcvcreux & Youna ( 1990) 
are all consistent to within 20%. In addition, the rest frame 160/^m 



in a single pixel, this detection seems marginal, but 
the local background of 0.03 counts/pixel implies 
that the likelihood of a random occurrence is only 
0.04%. The derived flux from BR1202N is fo.5 -2keV = 
4.8;^ jxl0-i^erg cm'^s^^ (95% confidence; iReeened 
implying Lo.5-2keV '-3xl0''*erg s'^. 

If the X-ray emission from BR1202N is real, its high 
X-ray luminosity suggests that BR1202N also hosts a 
luminous AGN. In addition, our new [C II] line detec- 
tion, large dust emission, and the presence of abundant 
molecular gas suggest that BR1202N is forming stars at 
a high rate. BR1202S and BR1202N thus form a system 
of two massive galaxies undergoing an interaction or a 
merger, separated by a projected distance of ~25 kpc, 
observed when the universe was only 1.2 Gyrs old. The 
presence of a highly obscured AGN (and an optically 
bright QSO) during a massive merger is consistent with 
the scenario predicted in recent galaxy collision simula- 
tions thMmvesti^a^^^ evolution of a central black hole 
(e.g. iHopkins et ^112001 . 

A high degree of clustering near massive galaxies is 
expected in the the oretical models of lar ge scale struc- 
ture formation (e.g. iWhite fc ReesI 1 197*81) . This is sup- 
ported by observations such as the higher surface density 
of submillimet er sources in the fiel ds surrounding high-z 
radio gal axies iStevens et al.l )2003*). Although an earlier 
study bv lGiallongo'eriLl.l l| 1998) found little evidence for 
clustering of galaxies near BR1202-0725 with R' < 25, 
two massive galaxies in proximity of each other makes 
BR1202-0725 a particularly interesting test ground for 
the hierarchical scenarios, soon after the reionization. 

5. SUMMARY 

We present a 3" resolution 900/im continuum image 
and a detection of the redshifted [C II] line emission from 
the z = 4.7 QSO BR 1202-0725 obtained using the SMA. 
This is one of the first detections of the [C II] line from 
a high redshift source. The [C II] line is associated with 
BR1202N. The low [C II]-FIR ratio of - 3.8 x lO"* is sim- 
ilar to local ULIRGs. X-ray emission is clearly detected 
from BR1202S, and at 99.6% confidence from BR1202N, 
suggesting that BR1202-0725 is the first example of a 
pair of AGN hosts at z ~ 4.7. 

The authors thank the anonymous referee for valuable 
comments, which have improved this work significantly. 

emission for BR1202-0725 traces the same warm dust (and hence 
the same M^j) as measured by IRAS I Yun et al. 2006). 
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Fig. 1. — {Left) 900/im continuum map of BR1202-0725 made by averaging line free channels in the LSB. The crosses mark the locations 
of the CO (2-1) peaks from Cari lli et af] 120021) . The contour levels are -2,3,4,5,6, 7,8,9,10,llt7 {la = 3.0 mjy beam-^). The beam is shown 
in the lower left corner. The SMA line profile of BR1202N {middle) and BR,1202S {right) obtained at the peak pixels using 120 km s~^ 
velocity averaging. Velocity of km s~^ corresponds to 334 GHz. The solid lines arc schematic representations of the central velocities 
and linewidts of the CO (5-4) line observed in BR1202N and BR1202S bv .Omont et al.. 1.1996i) . where the peaks have been scaled to match 
the [C II] line intensity in BR1202N. 
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Fig. 2.— The radio to FIR SED of BR1202-0725. Two values are given when BR1202N and BR1202S are spatially resolved. The soUd 
dark line represents the fit to the SED of the entire BR1202-0725 system, while the SEDs are not unique fits for BR1202N and BR1202S. 
The parameters used for the SEDs are T^j = 80K, a = 0.75, /3 = 1.75, f„th = 6.0 and SFR = 2000 Mq yr~^ for BR1202N and T^ = 65K, 
a = 0.75, /3 = 1.50, fnth = 2.0 and SFR = 4500 Mq yr~^ for BR1202S, where a, /3, and f„th are th e synchrotron spectr al index, dust 
emissivity index, and the normalization factor to account for the non-thermal synchrotron emission fsee lYun fc CarilH200l for details). 
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TABLE 1 
Summary of SMA Observations 



Date 




f 2 
('Source 


rp 3 

sys 


RMS phase'' 


02/22/2005 


5 


5.3 


240 


5 


04/05/2005 


6 


5.2 


310 


15 


03/21/2006 


7 


1.3 


300 


14 


04/19/2006 


8 


4.5 


250 


6 


05/11/2006 


5 


3.3 


300 


11 



"'"Number of SMA antennas used. 

"^On source integration time in hours. 

^The average DSB system temperature in Kelvins. 

^The average RMS phase (in degrees) taken over all baselines for each track. These are calculated toward 3C279 for the first two tracks and 3C273 
for the three tracks taken in 2006. 
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TABLE 2 
Summary of Luminosity Ratios 



Source 


LpiR (Lq) 


L[cii]/Lfir 


Rcf. 


Normal Galaxies 


10^ - 10" 


0.1 - 1% 


1 


ULIRGs 


~ 1012 


0.01 - 0.1% 


2 


High-z Sources 








BR1202N {z = 4.7) 


1.2 X 10^3 


0.04% 


3 


BR1202S (z = 4.7) 


2.6 X 10^^ 


< 0.03% 


3 


SDSS J1148+5251 (z = 6.4) 


1.2 X 10^3 


0.02% 


4 


CL 1358+62 (z = 4.9) 


2.4 X 10i2 


< 0.4% 


5 



References. — fll lMalhotra et alj 1200111 : (2l lLuhman et al.l r2003D : (3) this work; (4") IMaiolino et al.lt200S : fSi lMarsden et al.l f2005ll . 



